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In this study, we have investigated ground state properties of one- dimensional antiferromagnetic 
spin-1 chain with single-ion anisotropy at very low temperatures using the Transfer Matrix method. 
Magnetic plateaus, phase diagram, specific heat, susceptibility of the spin chain have been eval- 
uated numerically from the free energy. Results are in good agreement with the experimental 
data for the spin-1 compounds [Ni2(Medpt)2(/i-ox)(H20)2](C104)22H20, [Ni2(Medpt)2(At-ox)(^t- 
N3)](ClO4)0.5H2O, Ni(C2H8N2)Ni(CN)4 and Ni(CioH8N2)2Ni(CN)4.H20. However, spin-Peierls 
transition have not been observed in the temperature dependence of specific heat and magnetic 
susceptibility. 
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The physics of low-dimensional i.e., one-dimensional 
(ID) or quasi-one-dimensional (QID) spin-S (S> 1) 
chains has attracted a considerable amount of attention 
after the prediction by Haldane^ that a ID Heisenberg 
antiferromagnet should have an energy gap between the 
singlet ground state and the first excited triplet states in 
the case of an integer spin quantum number, while the 
energy levels are gapless in the case of a half integer spin 
quantum number. The most fascinating characteristic of 
these low-dimensional systems is that they show mag- 
netic plateaus i.e., quantization of magnetization at low 
temperatures near the ground state. This phenomenon 
has been observed not only in Haldane spin systems but 
also in other spin gapped systems for instance; spin- 
Peierls and spin ladders systems. A general condition of 
quantization of the magnetization was derived from the 
Lieb-Schultz-Mattis theorem^, for low-dimensional mag- 
netic systems. The fact that, Oshikawa, Yamanaka and 
AfHcck (OYA)^ discussed this plateau problem and de- 
rived a condition p{S — m) = integer, necessary for the 
appearance of the plateau in the magnetization curve of 
one-dimensional spin system, where S is the magnitude of 
spin, m is the magnetization per site and p is the spatial 
period of the ground state, respectively. The 25* -I- 1 mag- 
netization plateaus (contained the saturated magnetiza- 
tion m — S) can appear when the magnetic field increases 
from zero to its saturation value h^. Theoretical studies 
have suggested the realization of quantization of magne- 
tization in various systemS )^i'*i^i^i'''i^i^i^°i-^-^i^^i-^^i^'^ and it 
has been observed in experimental studies >i^*iSiiLi2iiS 

Recent years, antiferromagnetic (AF) spin-1 systems 
among the other low-dimensional gapped spin-S systems 
have drawn attention from both theorists and exper- 
imentalists in literature. So far, many QID gapped 
AF spin-1 systems, which are called Haldane systems, 
spin-Peierls and ladder compounds, were synthesized 
and to understand the ground state behaviors of the 
spin-1 AF Heisenberg chain have been extensively 
studied. The first QID spin-1 Haldane gap compound 
was Ni(C2H8N2)2N02(C104) (abbreviated NENP) 



which was synthesized by Meyer et aliSS in 1981, just 
before the Haldane's prediction appeared. Renard et 
al.^ showed that NENP had the energy spectrum 
predicted by Haldane; the magnetic susceptibility 
decreased steeply at low temperatures for all of the 
crystal axes suggesting a single ground state and 
neutron inelastic scattering showed that the magnetic 
excitation had a gap at the AF zone center. Other 
QID spin-1 compounds Ni(C2H8N2)2Ni(CN)4 (ab- 
breviated NENC))^^i^3,24,25 Ni(CiiHioN20)2Ni(CN)4 
(abbreviated NDPK),^^*^^ Ni(CioH8N2)2Ni(CN)4.H20 
(abbreviated NBYC),S2*2^ Ni(C5Hi4N2)2N3(PF6) (ab- 
breviated NDMAP) ,^6,27,28 g^^^ Ni(C5Hi4N2)2N3(C104) 

(abbreviated NDMAZ)— have also been identified as 
effective Heisenberg chains. This class of compounds 
exhibits a non-degenerate ground state which can be 
separated from the lowest excitation. As predicted 
by Haldane, it is experimentally shown that these 
systems have an energy gap between the single ground 
state and first excited triplet.^^'^^i^'^i'^'^i^^i^^i^^i^'^i^^ In 
addition, exception Haldane systems, several QID AF 
spin-1 gapped compounds such as AF dimer compound 
[Ni2(Medpt)2(/i-ox)(H20)2](C104)22H20,^ and AF 
alternating chain compound [Ni2(Medpt)2(/i-ox)(/^- 
N3)](ClO4)0.5H2O,l^ where (Medpt=methyl-bis(3- 
aminopropyl) amine), the ladder system 3,3',5,5'-tetrakis 
(N-tert-butylaminxyl)biphenyl (BIP-TENO),^'^i30,3i ^j^g 
spin-Peierls systems22iS*^ii2^ (Li,Na)V(Si,Ge)206 have 
been synthesized and magnetic properties have been 
studied. 

The magnetic quantization behavior has been ex- 
perimentally observed in only several gapped spin-1 
AF materials as predicted by OYAi^ For example, 
Narumi et alji^J^ observed a magnetic plateaus in 
the magnetization curve for both [Ni2(Medpt)2(/i- 
ox)(H20)2](C104)22H20 and [Ni2(Medpt)2(M-ox)(Ai- 
N3)](ClO4)0.5H2O. Goto et al}'' reported the existence 
of the magnetization a plateau at 0.25 in spin-1 3, 3', 5,5'- 
tetrakis (N-tert-butylaminxyl)biphenyl (BIP-TENO). 
On the other hand, in theoretical studies, Chen et 
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air- employed the classical Monte Carlo method to 
investigate the magnetization plateaus of ID classic 
spin-1 AF Ising chain with a single-ion anisotropy under 
the external field at low temperatures, and they showed 
that the system has 25' -I- 1 magnetic plateaus. Tonegawa 
et aliS observed the plateau of at m =0.0, 0.5, and 1.0 
in the ground state of spin-1 AF Heisenberg chain with 
bond alternating and uniaxial single-ion anisotropy. 
Sato and Kindo obtained a magnetic plateau at m =0.0, 
0.5, and 1.0 for spin-1 bond alternating Heisenberg chain 
using renormalization group methodi^^ 

Other thermodynamical features of low-dimensional 
AF spin-1 systems have been extensively studied in the- 
oretical concepts as well as in experimental studies. For 
example. White and Huse^fi calculated the ground state 
energy and Haldane gap of a spin-1 AF Heisenberg chain 
using the density matrix renormalization group tech- 
niques. Yamamoto and Miyashita'^^ studied the specific 
heat and magnetic susceptibility of a spin-1 AF chain by 
Monte Carlo simulation. The spin-wave theory was used 
by RezendeS to study the temperature dependence of 
excitation spectra, and it has been found that the exci- 
tation energy or Haldane gap increases with increasing 
temperature in the region of low temperature. Cam- 
pana et ali^ discussed the specific heat and magnetic 
susceptibility of a finite chain by a numerical calculation 
method. Bac'*'' employed a second-order Green function 
method to discuss the variations of internal energy and 
specific heat with temperature for a spin-1 Heisenberg AF 
chain. Li and Zhu,^^ using a self-consistent mean-field 
approximation, studied the thermodynamic properties of 
a spin-1 AF Heisenberg chain with an anisotropy in the 
presence of an external magnetic field at a finite temper- 
ature, they found that the internal energy and Haldane 
gap increase with increasing temperature. Batchelor et 
al4^ investigated the ground state and thermodynamic 
properties of spin-1 via the integrable su(3) model. Lou 
et aljH calculated the magnetic susceptibility and spe- 
cific heat of one-dimensional spin-1 bilinear-biquadratic 



Heisenberg model using transfer matrix renormalization 
group. On the other hand, a field-theoretic approach to 
the Heisenberg spin-1 chain with single- ion anisotropy 
was suggested by Tsvelik»M 

Now, it understood from experimental and theoretical 
studies that quantization of magnetization appears due 
to gap mechanism originating from dimerization, frustra- 
tion, single-ion anisotropy, periodic field so on. However, 
it notes that low-dimensional gapped spin systems have 
still open problems and noteworthy physics. Therefore, 
in this study, we interested in ground state properties of 
ID AF spin-1 chain. 

In general, the most of the theoretical studies about 
low-dimensional spin-S systems have been based on 
Heisenberg Hamiltonian because of non-trivial quantum 
effects. However, in several theoretical studies have been 
shown that it is possible to used classical spin systems 
to obtain magnetic plateaus and other thermodynami- 
cal behaviorsAiiii^ Considering this fact, we have con- 
structed herein the Hamiltonian of the ID spin-1 sys- 
tem in a quasi-classical manner which allows to use the 
Ising type variable instead of the quantum spin opera- 
tors. Such a Hamiltonian with single-ion anisotropy for 
ID spin-1 systems is described by 
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where J denotes the exchange coupling of antifcrromag- 
netic type {J > 0), D describe the single-ion anisotropy, 
and h is the external field. Also, S'f refers to spin of 
magnitude 1, which takes on 0, and ±1 values. 

We used Transfer Matrix method to obtain ground 
state properties of ID AF spin-1 system which was de- 
fined by Eq. (1). Using this method, it is possible to get 
the analytical and comparable results for ID infinite spin 
system, though it is very simple approach. 

For the Eq. (1) transfer matrix V is given as 
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where (3 — l/kT. Using relation (2), partition function 
of the ID spin-1 AF systems can be represented in terms 
of transfer matrix V as 



Z = TrV^. (3) 
Hence free energy per spin of this system is given by 



We are interested in the three thermodynamical expres- 
sions of (i) ferromagnetic order m, (ii) the specific heat 
C, and (iii) the magnetic susceptibility X; which respec- 
tively are 
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A conventional way to calculate free energy Eq. (4) which 
is required to obtain expression in Eq. (5) is that parti- 
tion function is expressed in terms of eigenvalues of trans- 
fer matrix (2) as 



(6) 



The eigenvalues of the Transfer Matrix (2) have been 
carried out with the aid of Mathematica 4.0 package. 
However, we have not listed herein since they occupied 
many pages. After we tested all eigenvalues numerically, 
we have sorted and so-called them from the biggest to 
smaller as Ai, A2, A3 respectively. Using standard as- 
sumptions of Transfer Matrix method, Eq. (6) was de- 
fined as, 

2 = TrF^ = Af[l + (^)^ + (^n. (7) 
Ai Ai 

It is clearly seen that on the RHS of Eq. (7) the second 
and the third terms goes to zero as ^ 00 since both 
terms are being | ^ | < 1 and | ^ | < 1. Consequently, 
the free energy Eq. (4) is reduced to 



f{h,T) = -fcTlnAi 



(8) 



then, some of the thermodynamical quantities of the sys- 
tem, we interested in Eq. (5), can be calculated numeri- 
cally using Eq. (8) which is depend on the biggest eigen- 
value Ai of the matrix (2). 

Obtained numerical results are followed: In Fig. 1(a) 
and (b), the magnetization m is plotted ai T — 0.01 
(units by J) as a function of external field h ioi D — 0.5, 
and D = 1.0, respectively. It is seen that 25* -I- 1 = 3 
plateaus appear (i.e. m =0, 0.5, and 1) for all values of 
positive single-ion anisotropy D at the low temperature 
near the ground state under external field h. No magne- 
tization appears in the low field region. By increasing the 
field there is a plateau of which starting point is called 
initial critical field he and then with increasing the field, a 
plateau of which starting point is called saturated field hs 
occurs again. These results clearly indicated that num- 
bers of plateaus of classical ID spin-1 AF system obey to 
OYA criterion as well as Monte Carlo prediction,'* and 
compatible with the experimental data for compound 

[Ni2 (Medpt)2 (m-Ox) (H20)2] (C104)22H20.15 

It suggested that single-ion anisotropy plays a signif- 
icant role for magnetic plateau.* In order to examine 
the effect of single-ion anisotropy on the magnetization 
plateaus, the magnetization m was calculated at finite h 
for a series of values of D (0.1 < D < 1.5), and the data 
were plotted in Fig. 2. For D > 0.0, three plateau lines 
placed at TO = 0.0, to = 0.5, and m — 1.0 are divided by 
the initial field and saturated field lines as shown in Fig. 
2. The longitudinal coordinate of the circular-dot line is 



the beginning point of the field for the appearance of the 
plateau to — 0.5 which corresponds to critical field he 
values, and its ending point is signed by the square-dot 
line which corresponds to the saturated field hs values. 
In addition, the distance between both the lines for the 
same value of D is the width of the plateau to = 0.5. 
For D > 1.0, the value of saturated field hs increases 
with the increase of D. The initial field he decreases 
with the increase of D in the interval 0.0 < D < 1.0 and 
increases for around D > 1.0. Almost same magnetiza- 
tion phase diagram for spin-1 has been carried out using 
Monte Carlo method with finite spin in Ref.4. 

Figure 3 (a) and (b) show the specific heat C plot- 
ted as a function of temperature T, for fixed value of 
single-ion anisotropy field D = 0.2 at various values 
of external field h — 0.3, 0.5, 0.7, and for fixed value 
oi h ^ 0.1 at various values of 13 = 0.1,0.3,0.5, re- 
spectively. These curves show a board maximum as 
Schottky-like round hump about between T = 8 and 
13. But, it is no means that phase transition occurs 
in one-dimensional system. Schottky-like round hump 
in the specific heat probably reflects energy of the sys- 
tem. For several higher external fields and single-ion 
anisotropy fields, at fixed value of D and h, respectively, 
the peaks get smaller and move towards zero tempera- 
ture. It seems clear that the peaks are related to the 
magnetization phase diagram of the ground state of AF 
Ising chain with single-ion anisotropy under finite mag- 
netic field. Furthermore, the specific heat curves of the 
spin chain which have obtained as a function of tem- 
perature are in good agreement with the experimental 
data for the compounds Ni(C2H8N2)Ni(CN)4t2ai2i212^ 
Ni(CioH8N2)2Ni(CN)4H20,2^iSi24 and the theoretical 
studies^iSLS^ in literature, respectively. 

An other way is to understand the ground state be- 
havior of the system is investigation field dependence of 
the susceptibility. In this reason, the susceptibility were 
plotted in Fig. 4 as a function of the magnetic field. It 
is seen that there are two peaks in the susceptibility for 
D — 0.5 at T = 0.01. First peak occur at initial field 
value, and second peak occur saturated field value for 
fixed value of D. Both of the peaks indicate the criti- 
cal fields where magnetization plateaus appear. Similar 
two peaks have been observed in experimental study for 
compound [Ni2 (Medpt) 2 (/i-ox) (H2 O) 2] (CIO4 ) 22H2 O 

On the other hand, for different values of the mag- 
netic field, the susceptibility was plotted as a function 
of the temperature T for fixed value oi D = 0.5 
in Fig. 5. The susceptibility curve for h — 1.0, 
marked by square-dot line, has relatively a sharp peak 
around T = 5. However, it has a broad maximum for 
h = 0.5, and h = 0.1 which are denoted by circular- 
dot and triangular-dot line, respectively. We roughly 
say that susceptibility curves decay exponentially 
with increasing temperature by the following relation 
x{T) w exp {—Eg/kBT). This characteristic behavior 
can be interpreted that AF spin-1 chain with the single- 
ion anisotropy has the Haldane gap like Heisenberg 
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AF systems. Our results clearly compatible with the 
experimental results for compounds [Ni2(Medpt)2(/i- 
ox)(H20)2](C104)22H20^5 [Ni2(Medpt)2(/i-ox)(/x- 
N3)](ClO4)0.5H2O,i^ Ni(C5Hi4N2)2N3(PF6),Si and 
3,3',5,5'-tetrakis (N-tert-butylaminxyl) biphenyljii and 
theoretical resultsSLSS in literature. However, in the 
present study, the specific heat and magnetic susceptibil- 
ity curves which obtained as a function of temperature 
have not a spin-Peierls transition. Because in the 
compounds synthesized experimentally, a spin-Peierls 
transition in C and x lo"^ temperatures are due to 
a small amount of magnetic impurities and derivations 
from stoichiometryi^ 

In the majority of plateau mechanism which have been 
proposed up to now the purely quantum phenomena play 
a curial role. The concepts of magnetic quasi-particles 
and the strong quantum fluctuations are regarded to be 
first important for understanding of these processes. Par- 
ticularly, for a number of systems it was shown that 
the plateaus are caused by the presence of the spin gap 
in the spectrum of magnetic excitations in the external 



field. However, we have used quasi-classical Hamilto- 
nian (1) in this study, and magnetic plateaus, magnetic 
phase diagram, specific heat, magnetic susceptibility of 
ID AF spin-1 chain have been obtained. Obtaining re- 
sults clearly consistent with experimental and theoretical 
results which ensure that quasi-classical approach can be 
used to examine the ground state properties of ID AF 
spin-S gapped systems as spin-1. 

Finally, we conclude that the single-ion anisotropy in 
the Ising model has a significant effect on the magnetic 
properties under the external field. Therefore, we studied 
this phenomenon in the one-dimensional spin-1 AF Ising 
chain with single-ion anisotropy. When the external field 
varied from the zero to the saturated field, the 2S-f 1 step- 
like plateaus occurred owing to the co-existence of AF in- 
teraction of {S^S^j^i) and positive single- ion anisotropy 
{Sfy. Hence, the energy levels, separated by the crys- 
tal field (i.e. single-ion anisotropy) in the system of two 
particles, correspond to the plateaus in the ground state. 

We acknowledge useful discussion with Hamza Polat. 
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FIG. 1: The magnetization m as a function of magnetic field 
h for (a) D = 0.5; and (b) D = 1.0. Where he and hs are the 
initial and the saturated field, respectively. (T = 0.01, unit 
by J) 




FIG. 2: Magnetization phase diagram of the ground state of 
antiferromagnetic Ising chain with single-ion anisotropy under 
finite magnetic field. The circle dot-line and the square-dot 
line represent the initial field and the saturated field, respec- 
tively. 
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FIG. 3: The specific heat C as a function of the temperature 
T (unit by J): (a) h = 0.3, 0.5, 0.7 for fixed value of L> = 0.2 
(b) D = 0.1, 0.3, 0.5 for fixed value of /i = 0.1, respectively. 




FIG. 4: Susceptibility as a function of magnetic field h for 
D = 0.5 and T = 0.01 (unit by J). 
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FIG. 5: Susceptibility as a function of temperature for various 
h and D = 0.5 (unit by J). 



